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I. On the Instability of the Pear-shaped Figure of Equilibrium of a 

Rotating Mass of Liquid. 

By J. H. Jeans, M.A., F.RS. 
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1. The main result obtained in this paper is sufficiently indicated by its title. In a 
recent paper # I showed that the stability of the pear-shaped figure could only be decided 
on after the figure itself had been calculated as far as terms involving the cube of the 
parameter 6, which is used to measure the deviation of the pear-shaped figure from 
the critical Jacobian ellipsoid. In the present paper these third-order terms are 
calculated, and the pear-shaped figure is definitely shown to be unstable. 

A statement of the results obtained, and a discussion of their bearing on the wider 
question of which this problem forms a part, will be found at the end of the paper 

2. The discussion has to begin with a determination of the potential of a distorted 
ellipsoid, carried as far as the third order of the small quantities involved. With a 
view to shortening very lengthy computations, it is convenient to arrange the algebraic 
solution in a form somewhat different from that previously given. The solution now 
given can readily be extended to any order of small quantities, and appears to lead to 
the most concise series of computations for terms of all degrees above the second. 

Potential of a Distorted Ellipsoid carried to the Terms of Third and 

Higher Orders. 

3. As before, the undisturbed ellipsoid (which will ultimately be supposed to be the 
critical Jacobian) is taken to be the surface A = in the family of surfaces 

J a 2 +\ 6 2 +A c 2 + A 
while the disturbed ellipsoid whose potential we require (which will ultimately be 
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supposed to be the pear-shaped figure) is taken to be the surface X = in the more 
general family, 



/+0 



x y z 

+ TT ~ + 



a 2 + X 6 2 + X c 3 + X 



+ 0-1 =0. ...... (2) 



Here is a function of x, y, 0, and X, different values for <j> representing different, 
distortions of the fundamental ellipsoid. 

If Y i9 V are the internal and external potentials at x, y, z (or their analytical 
continuations as explained in the previous paper, § 3), we have seen that V,-, V y must 
be of the form 

V<= ^(x)(/+0)dX, .......... (3) 

Jo 



where 



Vr(X) = 



V, 



yjs (x) (/+ d>) dX, (4) 

JK' ' 



irpabc 



(a 2 +x)(b* + X)(c 2 +X)} 1/ 



_ ^P 



abc 



A 



(5) 



In these integrations, x, y, z are treated as constants while X 7 is the value of X at 
the point x, y, z as determined by equation (2). Furthermore, cj> must not be selected 
at random ; equations (3) and (4) will only give the true potentials if <p is chosen so as 
to satisfy 



dX 3X/J 



xfr(x)v 2 <pdx+^(x) 



A dx " dX VBa;, 



0. . . (6) 



This equation, as before, is most conveniently solved by a solution 



<f> = u+fv, (7) 



in which u, v must satisfy the equations 



\^(X) V 2 (u+fv)dX + 4^(x)v = 0. 



• • « 



• . (oj 



4 (1+v) 



A 3^ 3\/ A A3a; 3x, 



,3fl3 das/ 



(9) 



4. We attack the second equation first. As in the previous paper, let us introduce 
new co-ordinates £ ?/, f defined by 

/y» i"y» 



a 2 +X A 
and the equation is found to reduce to (cf. equation (37) of the previous paper) 

«< i+ «>(£ + '!Hp(^'-* 



FIGURE OF EQUILIBRIUM OF A ROTATING MASS OF LIQUID. 

In this equation we may putjf+0 = 0, or, since <f> = u+fv> we may put 



l+v 9 



although of course it would not be legitimate to equate differential coefficients of these 
equal quantities. The equation reduces to 



o\ \l+v. 



u 



3£\l + y, 



V/J» • • • • « * V / 



and this may be readily solved in powers of the parameter e on assuming a series of 
the form 



u 



l+v 



= egi + e 2 g 2 + e% + e 4 # 4 + .... (11) 



Equating coefficients of the different powers of e we obtain 

dffi _ a 



^ = - iS -(2^-^j, and so on. 

To satisfy the first of these equations, g x must be a function of £ ^ £ only, say P. 
To satisfy the remaining equations, write 

"II Q A 1 

A = — — — , &c, so that 7p = 72 5 &c* 

a A dX A 

Then if P^ is written for 3P/3f, the equation for gr 2 is 

ax ~~ 4 ax r * ' 

of which the solution is 

<7 2 =-i(AP/+BP/+CP/)+Q, . . (12) 

where Q is a function of £ j?, f only. Proceeding in the same way, we find 



0s 



= |(A 2 P/P K + ... + 260?^^ + ...)-i(APA+ ...)+ R - • • ( 13 ) 



where R is another function pf £ j?, f only, and so on. 

B 2 
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5. We now return to equation (8). By a transformation given in the previous paper 
(p. 37), it is found that this is satisfied by taking 



v = when X = 0, . . . 



ciXltl 



VO f% WMIMlO t 



^ X dv dv 

A dx 3X 



= A — 
"*" 3x 



where a- may be any function of x 9 y, % and X, which vanishes when X 
X = X'. For equation (15) we may try provisionally a solution 



v 



w -\fw f +f 2 w" + . . . +f n w n + . . . . 



where w^iv\ iv" . . . are quantities satisfying 

\ A dx ' 3X 



-V 2 u-A 



ax 



• (14) 



> l x *J J 



and when 



(16) 



(17) 



A ex c\ 



O V vU* •••• • • • • • e • • i lO I 



At „ x dw n , dw. 



A dx 3X 



n+1 



V 2 w n _ u &c (19) 



After a good deal of simplification, the left-hand member of equation (15) is found to 
reduce to 



1 A3j; 3x, 



= -A 



'"N 



cX 



+ f- (/w' + 2/ V + . . . + V'X + • ■ • ) 



. I ZU I 



The quantity is so far undetermined. Let it be given by 



4 / u 



A\l+u 



(w f + 2fiv" 4- . . . + nf n ~hv n + . . . ) 



5 » 9 ♦ * » 



(21) 



so that the quantity in the square bracket on the right of equation (20) becomes 



~(f+ T ^-)(tv , + 2fw"+...+nf n - 1 w n +...) 



• « * « 6 



» \AA) 



When X — X' this vanishes through the factor f-\ . It will vanish when X = 

s 1+v 

through the last factor if we make 



w f — w n = w nf — ... = 0, when X = 0. 



(23) 



If this last condition is satisfied, expression (22) satisfies completely the conditions 
which have to be satisfied by a- in equation (15). Hence the value of v given by 
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equation (16) will be a solution of equation (15). Moreover, if w\ w'\' w f " ... all 
vanish when A = 0, v will also vanish when X = 0, so that equation (14) will also be 
satisfied. It follows that equations (16) to (19), with (23), contain a complete 
solution of equation (8). 

6. These equations can be solved in powers of the parameter e. Let us assume for 
u, v expansions in the form 

u = eui + £ 2 ^2 + e s u g + e% 4 + . . . , 
v = evx + e 2 v 2 + # 3 % + e%4 + . . . , 

and for w 9 w f jii) f \ ... &c, expansions 



4 4/1 



w = ei^i 4- e w 2 + e ^3 + 6 w i + • • * ? 
?i/ = eic/j + e 2 w' 2 + eV 3 + eV 4 + . . . , 



The coefficients in the expansions of u, v are of course not independent of those 
in the expansion (ll) already assumed for uf(l+v). We find easily enough the 

relations 

u i = 0i == -*• ? 

% = 03 + ^102 + ^201 ? 

^4 = 04 + ^103 + ^202 + ^01* &C. 

The value of 6 (equation (21)) is found to be 



& = 1 [efo (^ + 2fu/\ + 8/ V", + . . . ) 






+ e 3 {gr 2 (w , l + ^/w"i + 3/W" 1 + ... ) +& (w' a + 2/w'' a + ... )}]. 



On equating coefficients of different powers of e in equations (17)— (19), we obtain 



4 ^ A 3a; 3\ / 



— • V Wrjj » » • » * ♦ 



< 



V 



A dx 3x 



• « 



««- Y %1/^y CX>Q* } » * . * * 



(24) 



(25) 



r 



4 2 



'^^ + ^) = _W4A 



A dx 9\ 



8X 



^« + 2/w // 1 + ...) ', . . (26) 



< 



v. 



4(2f^ + ^- 2 )-4VH&c (27) 

"A ox dX ' 
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• 

There are an infinite number of such sets of equations, of which we shall need only 
the set for w%, W \, ...in addition to the two sets above. To simplify the equations as 
much as possible, let us limit ourselves to the type of distortion which leads to the 
pear-shaped series of figures of equilibrium. For this, as we saw in the previous 
paper, u x is of degree 3 in £ */, £ so that tv l must by equation (24) be of degree unity, 
and from equation (25) w\ must vanish. Similarly u 2 is of degree 4, so that w 2 is of 
degree 2, w\ is of degree zero, and w ft 2 = 0. Again u z will be of degree 5, w z oi degree 
3, w' 3 of degree unity, w'\ = ; and so on. 

The set of equations for w& v/ s . . . now reduces to 



r 



1 



V. 



\ A ox o\ / 
yj s " = w 3 w = . . . = 0. 



7. Let us now introduce the operator D, already used in the previous paper (§ 14), 
defined by 

D^A^ + Bp+CrTg . . (30) 

By differentiation with respect to X, we have 

3D i a 2 , i a 2 i a 2 



ax ~ a 2 df ' B 2 a* 2 c 2 af 2 

so that 3D/3X is the same as V 2 transformed into £ 17, £ co-ordinates. We can now 
solve equations (24) to (29) at once by transforming into £ */, £ co-ordinates, 
Equation (24) becomes in £ ?/, £ co-ordinates 

, a^ 8D 3D -p 

ax ax * 3x 

and since P is a function of £ ^ f only, this has the integral 

Wl = -iDP. • (31) 

No constant of integration must be added, for D vanishes when X = and w x must 
also be made to vanish when X = (cf. equations (23)). Thus equation (31) gives 
the true value of w x and we have also seen above that w\ — w f \ = ... = 0. 

The value of v x is accordingly 



v, = %V, 



-iDP=-i(AP«+BP„ + CPtf> (32) 
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From 55 6 and 4, we now have as the value of u. 



•ii 



u 2 = g 2 +v l g 1 

= -i(AP/+BP, a +CP/) + Q-iP (AP K +BP„ + CP ff ) 

= -iDF+Q. (33) 

The value of w 2 can next be found from equation (26). The right-hand member 
reduces to — V 2 u 2 , and the equation, expressed in £ ay, f co-ordinates, becomes 



dw 2 _ 3D 

ax ~ ax 



" *d\ v ax^ J 

of which the solution is 

w a = ,fe D a P a - iDQ . . . (34) 

and similarly equation (27) leads to 

^=-r,^D 3 P+AD 8 Q. (35) 

These values for w 2 and w' 2 are identical with those obtained in the earlier paper, 
although obtained by a slightly different method. We now proceed to the third order 
terms. 

On substituting for u s its value as obtained in § 6, we find, in place of equation (28), 

which, after a good deal of simplification, 

3w'o 



V 2 (flr 8 + V&) + V 2 («, a P) +>' 2 V 2 P + 4P 



ax 

= V 2 (g 3 + ^ 2 ) + 2s|- 2 P^ + K+>' 2 )V 2 P-4P^. . . . . (36) 

Introducing the various values which have been obtained for g& g 2 , v ly w 2 and w f 2 , 
we find, after simplification, that 

V%, + ^ 2 ) + 2sAp^ +Wa V*P 

= p A* d £ + qk? |? + "rAB ^ + S AB jj2 + . . . 

ax ax ax ax 

+ t*A — + P A — +crA —+...+ a— + (37) 

dX dX dX dX 
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where 



P = 



3 2 



Ld£ 



_ /ip3\ 



q = 



5 = 



t^ 



16 



3 2 



a 6 



40 2 



Lara? 



p ° /JLP 2 \ 



(m - p 



40 2 



dg% 



P) 6 P, 6 

( -P 3 ) — P 



afVar 



awaf 



(*F) 



(*P 9 ) 



tX'O* 



4 



idi% 



(PQ) - P 



3< 



a^a,,- 



Q 



, &C, 



ot = R^, &c. 



Relation (37) can accordingly be put in the form 



3 rj 



V 2 (9s+v ld2 ) +22ip^ + w 2 W = ^ [^D» (*F) - ID 2 (PQ) + DR + SP^'J. 



Equation (36), after transformation to £ jy, f co-ordinates and integration with 
respect to A, gives 



w 3 = -i[^D s aF)-iD»(PQ)+DE]-«/ 8 P-i .MVPdx.. . (38) 






Let us now put 



f w' 2 V 2 Pc£\ = ^A, (39) 

Jo 



so that A will be a function of A only. We then have 



*' <*' 3A 



'A' 



! 



>' 2 V 3 P^X = | fi~d\ =fiA-€(f+n a +?)l*dK 



. . (40) 



whence, on operating with V 2 , 



V6 



r\' 



fu/ a v 2 p d\ = i 



6 _2 2 

3A HA B c 



rV 



AcZA. 



* • » 9 I Til 







Equation (29), transformed to £ */, f co-ordinates, now becomes 



3w/. 



3A 



3 __ 



ivy, 



_l_ 

3 2 



A ^ D» (*I») -| ^ D 2 (PQ). + ^ DP' 



ax 



+ "S" gT (£ ■ A ') + lh£ gT 



fi 2 2 
ABC 



•f- 



At2A 
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giving on integration, 



w's = h [ 3 k D* (iF) - A D 3 (PQ) + W^l 



+ i£A + -h£[j^ + | + §)£ A d\. 

8. The solution is now complete as far as third order terms, but can be expressed 
in a more convenient form. We have found for the whole value of v 3> 

v s = Wz+fw'z 

= -i [ifc D 8 UP) - iD 2 (PQ) + DR] -iv> 2 Y 

+ <kfh^(m-^(PQ)+±im] + Z, .... (42) 

where Z is formed of terms involving the function A, and has its value given by 

Jo 

The term Z in v ?t gives rise to a term e^fZ in 0, and this in turn leads to a term 

+ (x) eyz = - ^ e»/Z, 

in the function $ (see §§4, 11 of the previous paper), from which the whole solution is 
derived. Using the value for Z which has just been obtained, we readily find that, in 
x 9 y, z co-ordinates, 

_ZE"te#fZ= -iTrpabce' 4t J?t f A ^i- • .... (43) 
A . o\ LA A Jo J 

We found, however (see footnote to p. 32 of previous paper), that for a given 
potential problem, the value of <& is not unique. If any function <£ gives a solution of 
the potential problem, then it was found that any other function of the form 

<I> + JL {f [\d\\, (44) 

9x L Jo J 

will give a solution of the same problem, provided that F is any function of x 9 y 9 z 9 and 
A, which vanishes when A = A' (i.e., when x, y 9 z 9 and A are connected by relation (2)), 
and u is any function of A whatever. Consistently with these conditions we may 
take 

<w = A, F = ^Trpabce 3 * ^ *t > 

and the new solution (44) becomes 

* +Wa ^p^i£j; AA } ( 45) 

VOL. GCXVII. — A. C 
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In this solution, $ already contains a term involving A, namely that given by 
expression (43). Combining this with the remaining term in expression (45), we find 
that the new solution can be put in the form 



u Teh T* 



r\' 



A d\ 



(46) 



where $ is the old solution <£ with the terms in A omitted. The last term in expres- 
sion (45), being proportional to e n (j>, is of the fourth order of small quantities. Thus 
in a solution as far as e s only, this term may be omitted, and $ = $ will be a solution. 
In other words the term Z may be omitted entirely from equation (42), and the 
remaining terms will still give an accurate solution for v 3 . 
Omitting this Z-term, we obtain for the third order terms, 

ih+f»s = Us+f(w B +fv/ z ) 

= AD f -(iP)-iD(PQ)+R 

- ifi nVD 3 (*F) - 1 D 2 (PQ) + DR} 

+ A/MTrkD 4 (*P 8 )-AD 8 (PQ) + *DTl} (47) 

This completes the solution of the general potential problem. 



Potential of the Pear-shaped Figure, 

9. Collecting the results obtained in §§ 3-8, we have found that as far as terms of 
the third order of small quantities, a value of <p which satisfies the necessary 
differential equation (6) is 

= e (u x +j%) + e 2 (u 2 +Jv 2 ) + e 3 (u s +fv 3 ) 
= e[P-i/DP] 

+ e 2 [Q-|DP 2 +/{^D 3 (P 2 )-~iDQ} +/ 2 { -xsWD^+^D^}] 
+ c 8 [^-i/DPQ + T feD a P-i/{^D 8 P-iD 2 PQ + DR} 



(48) 



At the boundary A = 0, this value of $ reduces to 



= eP + e 2 Q + e 3 R, (49) 

and since P 5 Q, R are entirely at our disposal, this is capable of representing the most 
general displacement possible, as far as the third order of small quantities. We have, 
however, to save the printing of additional terms, already assumed that P is of a 
degree not higher than the third in x, y> z and Q of a degree not higher than the 
fourth. Subject to these limitations, the potential of the ellipsoid deformed in any 
way can be obtained by inserting the value (48) for <j> into equations (3) and (4). 
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10. We shall now limit ourselves to the distortion which was found, in the earlier 
paper, to lead to the pear-shaped figure of equilibrium; For this particular distortion 
we found (cf. equations (93) and (94) of the earlier paper), 

P = ^(a^ + /3, 2 + yf s + ic),. . (50) 

Q = i[L^ + M^+Nf 4 +2^ s + 2m^ S +2w#V+2 (pf + q^+H 2 ) + s]. . (51) 

We shall now assume for R the value 

R = i-£[&t+Mn i +m i +2hV+2mi 2 f+2n£\ 2 +2 (Vf +&+&)+*], (52) 

this being found adequate to satisfy all the conditions for a figure of equilibrium of 
a rotating liquid mass. 

11. The values of u 1 +fv l and of u 2 +fv a were calculated in the former paper (see 
equations (121), (122) and (132)). In the same way the value of u g +fv 3 as given by 
equation (47) is found to be 

i h + f v . 6 = f [1 5fAV + 2f AB a 2 /3 + 2-f AC« 2 y + f B 2 a /3 2 + iBCa/3y + f C 2 ay 2 ] 

+ fV [l3|AV/3 +7iAB a 2 /3+2iAC« / 8 y +l|B 2 ^ +|BC/3 2 y + |C 2 /3y 2 ] 

+ fi 2 [l3|AV y + 2|ABa/3y + 7p.C a y 2 + fB 2 /3 2 y+fBC/3y 2 +l|Cy] 

+ t V [ l|A 2 a /3 2 + l|AB/3 3 + fAC/3 2 y] 

+ £? [ l|A 2 « y 2 + fAB/3y 2 +l|ACy s ] 

+ inT [ 3fA 2 a /3y + 2iAB/3 2 y + 2iAC/3 y 2 ] 

+ ?k [l3iAV + 2fAB a /3 +2|ACay + |B 2 /3 2 +iBC/3y + fC 2 y 2 ] 

+ & 2 k [ 3fA 2 a /3 +2iAB/3 2 + |AC/3y] 

+ g K [ 3fA 2 ay + |AB/3y +2iAC r 2 ] 

+ £*? [ l|A 2 a + |AB/3 + fACy] 

-f [ 2|ALa + £B(Lj8 + 2wa) + iC(Ly + 2ma)] 

-£V [ ljA(L/3 + 2wa) + fB(Ma + 2»/3) + JC (Za + m/3 + fiy)] 

-ff [ liA(Ly + 2ma) + iB(la + ml3+ny)+ |C(Na + 2ny)] 

~£n [ |A(Ma + 2w/3) +1|BMj8 + iC(M y +2Z/3)] 

-#* [ fA(N« + 2n y ) + iB(N/3+2Zy) +l|CNy] 

- W [ |A («« + m/3 + ray) + |B (My + 2Z/3) + fC (N/3 + 2Zy)] 

-f 3 [ liA(I*+2ap) • + iB(n K +aq + /3p)+ £C (m* + y£> + ar)] 

-^ 2 [ |A(^ + ag+/3p)+ |B(M/c + 2/3g) + £C .(k + y g + /3r)] 

-ff 2 [ |A (m/c + yp + ar) + £B (L + y q + /3r) + |C(N/c + 2yr)] 

-i C |A(2^+s«) + iB(2g« + Sy 8) + |C (2r* + sy)] 

+ i£ C&f 4 + Wn + W + 2 V£ 2 + 2tttf »# 2 + 2n£ V + 2 ( JJf + q>? 2 + rf ) + * 

2 
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- ff [2 6iA«« 8 + AB 3 /3 3 + r \Cy + 6 ^A 2 B a 2 /3 + 6 ^A 2 Ca 2 y 

+ l|AB 2 ay 8 2 + l|AC 2 ay 2 + liABCa/^y+x^C/SV + ABC^y 2 ] 

- & s / [6M 3 « 3 ^ + 5f A s Ba/3 2 + 1 f A 2 Ca/3y + 2 |f AB 2 /3 3 + AAC 2 /3y 2 + l|ABC/3 2 y 

- ^[6 A-AVy + if A 2 Ba/3y + 5|A 2 C«y 2 + tVAB 2 /3 2 y + 2 j-f AC 2 y 3 + l|ABC/3y 2 
-if [6AAV/C + lfA 2 Ba/3,c + lfA 2 Cay,c + T %AB 2 /3 2 K + -&AC 2 7 2 K + fABC/3y,c 
+ ff [6x 9 BA 2 L a + AB 2 (Ma + 2»/3) + AC 2 (Na + 2n 7 ) + |AB (L/3 + 2na) 

+ fAC (Ly + 2 ma) + |BC (la + m/3 + ny)'_ 
+ inf [if A 2 (L/3 + 2»a) + 2if B 2 M/3 + r 3 B -C 2 (N/3 + 2/y ) + l|-AB (Ma + 2«/3) 

+ |AC (Za+m/8+«y) + fBC (My + 2Z/3)] 
+ ^ 2 /[i|A 2 (Ly+2ma)+ t \B 2 (My + 2^) + 2t|C 2 Ny + fAB (la + mfi + ny) 

+ 1|AC (Na + 2ny) + fBC (N/3 + 2/y)] 
+ if [«A 2 (L, + 2ap)+ AB 2 (M, + 2 / 8g) + AC 2 (N K +2yr) + fAB(n, + ag + ^) 

+ f AG (niK +yp + a.r) + gBC {l/c + yq + /3r)] 
-f/[li!LA+ InB + imC] 

-tff [fltA+fJlB + 11C] 
-&'f [fntA+ llB+fNC] 
-if [fpA+ |qB+ £tC]+£fG, 



• (53) 



where/stands for Af 2 + Ly + C£ 2 — 1, and G is given by 



G = A {iyJ>AV+ \% 1 A 2 Ba 2 j 8+ -V°/A 2 Ca 2 y + ff AB 2 a/3 2 

+ i i AC 2 ay 2 + MABCa^y + ^B'/S 8 + if C V 
+ ^B 2 C/3 2 y+ABC 2 /3y 2 } 

- -^A 3 La - if A 2 B (L/3 + 2na) - if A 2 C (Ly + 2ma) 

- B %AB 2 (Ma + 2w/3) - e^AC 2 (Na + 2w y ) - 3 %ABC (la + m/3 + n y ) 

-ifB 3 M/3-MC 3 Ny- 6 ^B 2 C (My+2^/3)-^BC 2 (N/3+2/y) 

+ MA 2 2,+ s WJ« + AC 2 ^+TVBCl+ 1 %ACm + AABtl. . . (54) 

12. The potential of the distorted ellipsoid 

/-heP + e 2 Q + e 3 E = 0, ......... (55) 

can now be written down as far as terms in e 3 . The terms in e and e 2 , which of 
course involve only P and Q, have been calculated in the previous paper. The terms 
in e 3 in V 6 , the potential at the boundary X = 5 are (cf. equation (3)), 



e 



"00 

y/r (x) (u H +fv 8 ) d\ = —irpabc e s 

o 



/•oo 







u 3 -\-fv< 






*(£\, 



, . (56) 



in which u$+fv d has the value just given in equation (53), 
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Let us put 



■GO 







u **f v * d\ = x(t u x" + t 22 yU~t 33 z k + t 12 xY-{-t dl z 2 x 2 +t 23 yV (57) 



A 



The values of the various coefficients in this expression are found to be as follows :- 



'qo 



(* — 

VII — 



1 



'11 



AA 



[1 5fAV + 2f ABa 2 /3 + 2f ACa 2 y + |B 2 a/3 2 +iBCa/3y 

+|C 2 ay 2 -2f ALa-fB (L/3 + 2raa)-fC (Ly + 2ma) + iiL] dX 



>oo 







J^[26iAV + T W/3 3 + AC y + 6 ^A 2 Ba 2 /3 + 6 AAC« 2 y 

+ l|AB 2 a/3 2 + if AC 2 ay 2 + l|ABCa/3y + fVB 2 C/3 2 y + f 6 BC 2 /3y 2 ] cfo 



*oc 



+ I ~[6 T 9 B A 2 L a +AB 2 (M a +2^) + AG 2 (N a +2ny) 



+ f AB (L/3 + 2na) + f AC (Ly + 2ma) + |-BC (la + m/3 + wy)] dX 



'CO 



^j[iiAH+|Bn+iCm] + £^ dx. 



• • 



. . . (58) 



vi o — 



1 

rrroa [ ! 3H V0 + 7|ABa/3 2 + 2p.Ca/3y + If B 2 /3 3 + |BC/3 2 y 
o AA d 

+ |C 2 /3y 2 - l|A (L/8 + 2«a) - f B (Ma + 2n/3) -|C (Za + m/3 + ny) + |tl] dX 



»oo 



1 



AA 3 B 



= [26iAV + AB 3 /3 3 + ACV + 6 &A 2 Ba 2 /3 + 6l \A 2 Ca 2 y 



+ if AB 2 a/3 2 + if AC 2 ay 2 + liABC«/3y + f^y + ^BC^y 2 ] dX 



•CO 



1 



AA 2 B 



2td2 [6 A A3 * a /8 + 5f A 2 Ba/3 2 + if A 2 Ca/3y + 2LfAB 2 /3 3 



AAC 3 /3y 8 + l|ABC/3 2 y] dX 






CO 



AA 2 B 



[If A 2 (L/3 + 2na) + 2if B 2 M/3 + -ft-C 2 (N/3 + 2Zy) 

+ l|AB (Ma + 2n@) + |AC (Za + wi/3 + fty) + f BC (My + 210)] dX 



•CO 



+■ \ ^g[6AALa + AB 2 (M« + 2^) + AC 2 (Na + 2 W y) 



+ |AB(L/3 + 2na)+-fAC (Ly+2ma) + fBC (la. + m^ + ny)'] dX 



(oo 



1 



AA 2 B 



2 



[fAtt+fBm+|Cl]c£\ 



poo -i 

Jo ^A 3 BC 1 i A1L +i B «+i Cm 3^+ 2 



•oo 



G 



AA 2 B 



dX. 



• • • • 



. . . (59) 
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The value of C 13 can be written down from symmetry. 



f*CO 



'22 



AAB 



, [if A 2 a/3 2 + if AB/3 2 + f AC/3 2 y-fA (Ma + 2«/3) 



lfBM/3-|C (M y + 2^) + iJ«] d\ 



'CO 



1 



ZaxxjL) 1 



[6 AAV/5 + 5f A 2 Ba/3 2 + If A 2 Ca/3y + 2ifAB 2 /3 3 



+ AAC^y 2 + l|ABC/3 2 y] dX 



+ 



'QO -I 

o aTb 3 [it a " (L/3 + 2na) + 2 ^ WM P + A° 3 W + %) 

+ 1|AB (Ma + 2/1/8.) + |AC (Ja + m/8 + «y) + |BC (My + 2Z/3)] cfo 



foe 



1 



AAB : 



[|Att + fBJl + iCl] + 



f*00 



Gr 



AAB 



- C& A ■ • • » V ^ ^ / 



The value of C 33 can be written down from symmetry. 



^23 ~ 



-r-r%^ 2 [3fA 2 «/3y + 2|AB/3 2 y + 2iAC / 8y 2 
o AA.1J O 



|A(/a + W/3 + ny)-tB(My + 2Z/3)-fC(N/8 + 2^y)+|l]dX 



2/93 



-—L^ [6 AAV/3 + 5f ABa/3 2 + if A 2 Ca/3y + 2H AB 2 /3 



+ AAC 2 /3y 2 + |ABG/3 2 y] d\ 



»oo 



1 







AABC 



12 



[6-MVy 4- If A 2 Ba/3y + 5|A 2 Gay 2 + AAB 2 /3 2 y 



+ 2fgAC 2 y 3 + |ABC/3y 2 ] dX 



/»00 






1 



AAB 2 C 



.00 



+ 



1 







AABC 2 



[i|A 2 ( L/3 + 2na) + 2f |B 2 M/3 + T %C 2 (N/3 + 2Zy) 

+ |AB(Ma + 2»/8) + fAC(7a + m J 8 + ny) + PC(My + 2Z/3)]dX 

|-|A 2 (Ly + 2ma) + ^B 2 (My + 2Z/3) + 2lgC 2 Ny 

+ |AB (Za + m/3 + »y) + fAC (Na + 2wy) + |BC (N/3+ 2Zy)] dX 



o ^^[fAm+iBl+fCN]^x 



,00 



1 



AAB 2 C 



[|An+|Bfft+iCl]dx 



+ 2 



G 

AABC 



<iX. 



» • » » 



• • * 



»»»»»« «? 



• ( 61 ) 
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The values of the four coefficients Xl u tf 2 , tf 3 , tf 4 are most conveniently expressed in 
the form :— 

G 






/•CO 



*co 



So = <L — 



to* = & 



o AAB 

1 



»oo 





«oo 



Zaxxvv 



(fpA + iqB + |rC) d!x + £q 
(fpA + iqB + |rC) dx.+ Jr 



• GO 



o AAx> 



>oo 



dX 
o AAC W 



/ji 



2 






G 
o AAB 



— dX. 





► 00 



** = ^ + TT (tP A + il B + £ rC ) dx + is 

Jo /AJ± 



• 00 



*J0 Zaxx 



/•00 



/•GO 



G 







AAC 

G 

AA 



dX, 



dX, 



dX, 



(62) 
(63) 
(64) 
(65) 



where S u S 29 S 39 4 are quantities which do not depend on p, q, t, JB, being given by the 
following equations : — - 



A 



C?9 



,00 ^ 

o AA 1 K 1 3|AV + 2|AB a /3 + 2|AC a y + tB^ + iBC/3 y + fCy) 

— l|A(L/c + 2ap) — |B(«/c+ag + /3p) — |-C (mic+py+ra.) 
. + 2 6|A 3 a 3 + ftB'/S 8 + ftC 3 y 3 + 6ftA 2 Ba 2 /3 + 6 ft A 2 Ca 2 y 
+ if AB 2 a/3 2 + l|AC 2 ay 2 + liABCa^y + AB^y+x^BC^y 2 ] dX 

[6ftAV + if A 2 B«/3 + If A 2 C«y + ftAB 8 /? + ftAC V + fABC/3y] dX 



•oo 


/.GO 



K 










AA' 



^ [6 &A 2 La + AB 2 (Ma + 2nPY+ T 3 ^C 2 (Na + 2my ) 

+ |AB(L/3 + 2wa)+|AC(Ly + 2ma)+iBC(Za + m/3 + wy)]d!x 

+ I "-^ ft IA 2 (Lc + 2*p) + T p 2 (M* + 2/3q) + AC 2 (N* + 2 y r) 

+ f AB (^ + olq + /3p) + fAC (mK + yjp + ar) + |BC (Ik + yg -f /3r)] c£\ 

-^[liAfi+iBn+iCm]^. (66) 

o AA 



+ 



! o 



.00 



• GO 



1 



o AAB 5 [ * (3f AV + 2iAB^ 2 +fAC^y) 



-I 



as 


00 



K 



AAB 
1 



-f-A(te + a g + /3p)-fB (M/c+2/3g)-|C (h + yq + pr) 
+ 6ftAV/3 + 5f A 2 Ba/3 2 + if A 2 Ca/3y 
+ 2i f AB 2 /3 3 + ft AC^y* + f ABC/3 2 y] dX 

[6 T 9 ff AV + if A 2 B a /3 + if A 2 C ai 8 + ft AB 2 /^ + ftAC 2 y 2 + -|ABC/3y] dX 



AAB 



,»oo 



+ 



+ 



1 
AAB 



' 1 
o AAB 5 



if A 2 (L/3 + 2wa) + 2if B 2 M i 8 + T %C 2 (N/3 + 2Zy) 
+ lfAB (Ma + 2W0) + fAC (Za + m/3 + n y ) + fBC (My + 21$)'] d\ 
[|f A 2 (Lk + 2ap) + ftff (Mk + 2yQg) + T 3 e C 2 (Nk + 2yr ) 

+ f AB (n/c + ag + ftp) + fAC (m/c + yp + ar) + fBC (Ik + yq + fir)] d\ 



[fAtt + fBJl+|CI]cZX. 



(67) 
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The value of S 3 can be written down from symmetry ; that of S 4 , which is of the 
same general type, will not be required in the present investigation. 

13. This completes the solution of the particular potential problem which we have 
had in hand. It might naturally be feared that some mistake might have been made 
either in principle or in detail, and it must be remembered that even one mistake 
might invalidate the answer to the whole problem. I have, therefore, both here and 
elsewhere, taken the utmost care to check the accuracy of my work in every way. 
The following will, I think, show that no error need be feared in the solution which 
has just been obtained. 

The value of Y i9 if obtained accurately, ought to satisfy V 2 V,- = -~47r/), and the term 
— 47rp must come entirely from the terms independent of e in Y { . Thus the terms in 
Y { which are multiplied by e, e 2 and e 3 ought separately to be spherical harmonics. 

It was verified in the previous paper that the terms in e and e 2 were in actual 
fact of this form. The terms in e 3 will be harmonic if the quantity on the right of 
equation (57) is harmonic, and the conditions for this are expressed by the equations 

10Cii + C 12 + Ci3 = 0, (68) 

®*22 "T" " vj 2 i ^23 == v/, . . . y\)i7 J 

^ V33 Hr O Cj3 "T V23 =: U 3 \l U J 

3&i+& a +lr 8 = 0. .......... (71) 

I have inserted the values just obtained for Cn, C 12? ... 6i, ... in these equations, and 
have verified that they are all satisfied. (The necessary transformations of the various 
integrals are tedious, but involve no special difficulties.) 

It follows that the solution we have obtained gives accurately the potential of some 
solid of uniform density p. By the method explained in the earlier sections of the 
previous paper, it is easy to work the problem backwards and to verify that the 
equation of the boundary of the solid in question is obtained by putting X = in the 
equation/ + ^> = 0. Thus Ave verify that our solution gives accurately the potential 
of the solid of boundary 

/+ eP + e 2 Q + e 3 B = 0. (72) 

Conditions that Pear-shaped Figure shall be one of Equilibrium for a Rotating 

Liquid. 

14. In order that the figure determined by equation (72) shall be one of equi- 
librium for a rotating liquid, the potential at the boundary plus |V (x 2 + y 2 ) must, as in 
equation (98) of the former paper, be identical with 

— TrpabcO {f + eP f e 2 Q + e 3 R } + a constant. 
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Let us limit ourselves to the terms in e\ Using the value of Y b given by equation 
(56) and that of R assumed in equation (52), we find that we must have 

- ir P abce?x (c u as 4 + t^if + C 33 z 4 + t r2 x V + C 31 z V + t 23 y 2 z 2 + tf x x 2 + %y 2 + %z 2 + fc 4 ) 



a 



a 



irpabcOePi ^ ( & ~ + Jl | + H ^ + 21 f^ 4 + 2m 



c 



.2^2 



6V 






+2n sS +2 (^ + ^ +t l) +8 



On equating coefficients we obtain 



C 



ii 



4 Tio^L ; 



l r 

^" 2 aW 1, 



c, 



22 



'31 



4 a 2 & 8 ^'' 


C33 — 


4 «V '" 


<9 

2 aV lj 


t 12 = 


1 n 

2 aW* 1, • 



e •» 



(73) 



• • (74) 



^1 — 2 rP ' 



a 



% 



1 



2 27 4 

a - 



q; 



^3 ~ 2"~2~4^ .... (75) 



4 a** 



(76) 



These equations, in addition to those found for the first- and second-order terms in 
the previous paper, express the condition that the third-order figure (72) shall be a 
possible figure of equilibrium. 

15. On substituting the values of C u , C 12 ... C 33 which have been obtained in § 12 into 
the six equations (73) and (74), we obtain a system of six equations from which it is 
possible to determine the six unknowns 5L, ftt, N, I, tit, U+ The solution is actually 
effected in § 17 below. 

If we substitute the three values of tr l5 % and tf 3 obtained in § 12 (equations 
(62)-(64)) into the three equations (75), we obtain three equations which can be 
written in the form : — • 



1 e_ , 3 

,6 



,00 



/»0C 



V[t a« + i hA&~ 2 



i^ +iq 



/»CC 



AA 2 ' 4 Jo AA' 



oc 



G 



5 = ^-2! ^- 2 dA, . . (77) 



»OD 



3 



JBlU/A 



e 



,o. 



dX 



/%CC 



1 



IP I TTr5" + qi2:^rz+4 j ^jg~2 



AAB 



VcfV 



dX 



AAB 2 , 



+ir 



QdX 
AAB 



/*GG 



4— 2 



G 



J-0 



AAB 



dX, (78) 



.oc 



fp 



/•OO 







AAC 



H 



dX 



e 



.00 



GdX 



.00 



AAC 



+ rlt ^ +4, oAAC _2 



dX \ 



»oo 



aV 



AAC7 



S,-2 



G 



AAC 



dX. . (79) 



It will be seen that p, £(, t do not occur on the right-hand sides of these 
equations. 
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These equations are not, as might at first be thought, a system of three simple 
equations determining p, q, r. They will be found to be of the type known as 
" porismatic " # ; that is to say, equations which are inconsistent unless the coefficients 
satisfy a certain relation, and such that, when this relation is satisfied, the equations 
have an infinite number of solutions. 

Let us, for brevity, write the equations in the form : — 



k$ + h\t{ + W X X = % 






>- 



6 * 



. . . . (80) 



j 



Let us use also the abbreviated notation of the previous paper (p. 50), such that 



.00 



Ci 



XdX 



Jo 



A ABC 



on > 



c 2 = 



XdX 

AA O 



Co = 



Xc?X 



Then, by simple transformations of the integrals, we obtain 



&i 



&„ = 



Aye 



1 


/20 

w 


Aa 2 


3 

4 a 3 


C 3? 


3 


C 2? 



— Co — c 



3 /? 



¥, 



1 



Ay 1 



20 



l_ 

4& 2 CiS ' 



F, 



1 



"A 



4c 



2 C 2) 



2 



46 2 \a 2 6 2 



3c 3 — Cj ), & 2 — ,, 2 c i» 



4a 



^ 



1 



4fr 



Ci, 



&" s 



1 / 20 



2 \ ^2^2 



4<r Wc 



4<T 



3co— c 



> . . . . (81) 



2 ^1 /j 



With these values of the coefficients, it will be found that equations (71)- (7 3) of 
the previous paper reduce to 



so that 



V 



2\a + 2k\fi + 2F 1 y = 
2A 2 a + 2&' 2 /3 + 2F 3 y - 
2Jc s ol + 2# a /3 + 2&" 3 y - 

= 2a, - q = 2/3, r 



» • e 



» » • 



(82) 



2y, 



(83) 



is a solution of equations (80), when % x — % 2 = % s = 0. 

Without this detailed inspection of the equations, it could have been foreseen that 
this would necessarily be the case. For our general solution 



= eP + e 2 Q + e 3 R, 



e • » P » 



. . . (84) 



* See ' Hobson's Plane Trigonometry,' § 73, or Wolsten holme, ' Proc. London Math, Soc.,' vol. 4. 
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must cover all possible figures of equilibrium as far as the third order. One such 
figure is, however, known to be 



= e.O + e 2 .O + e 8 P, . . . 



• 



e e • • I O «J / 



and the corresponding solution for JJ, q, X is that expressed by equations (83). 
Prom equations (82), it follows at once that we must have 



k u k 2 , k 1 

~k V V 1 

A/ 3 , th 3 , /I 3 



- 0. 



• • a • 



. . (86) 



Indeed, it is now clear that this is precisely the equation which determined the 
existence of the point of bifurcation on the series of Jacobian ellipsoids.* 

If this were the only relation between the coefficients in equations (80), -these 
equations could have no solution other than p = q = X = oo . Let us, however, multiply 
the three equations (80) by the three minors of k" l9 k f/ 2 > k f/ 3 in the determinant of 
equation (86), and add. We obtain a relation of the type, 



p.o + q.o + r.o 






A/ 3 j th 3, 



• (87) 



3 



and it is clear that these equations can now have a solution in which p, q, X are not all 
infinite if we have 



Jr V 



A>9, A/ 



'2? 



2? 



0. 



fCo. IC >: 



'3> 



3) 



(88) 



It is only when this relation is satisfied that it is possible to continue our linear 
series of equilibrium configurations beyond the second order terms. When it 
is not satisfied, our third order solution (84) lapses back into the solution (85), 
namely <f> = e 3 P, which is virtually the first order solution with e s replacing e as 
parameter. 

16. The relation (88) can be expressed in a much simpler form. 

Independently of the values of }J, q and r, we have already seen (equation (71)) 
that we must have 



3)*! + & 2 + & 8 = 0. .......... 

* As to the relation of this equation to the general theory, see § 36 of my previous paper. 



(-89) 
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On substituting the values for fr l5 tr 3 and % given by equations (62)-(64), ' and 
equating coefficients of p, C{, T in equation (89) to zero, we obtain 



3M 1 + % 2 + » 3 - 0. 



3*! + £ 2 + *« - § - 6 






3Fj + F 3 + &" 3 



-JL 1 

" 2 a V * 



With the help of these relations equation (88) reduces quite simply to 



\a b / \a I) 



^(^tf 2 - ll -k 2 ) = %(^kr i - 1 ^k 1 ). ...... (90) 



In this equation the coefficients of 3^ and % 2 depend only on a, &, c the semi-axes 
of the Jacobian ellipsoid, and so are fully known. The quantities % x , % 2 however 
depend on the second-order coefficients L, M, N, ... p, q> r, s. These were calculated 
in the previous paper, but p, g, r, 5 could not be fully determined, since they were 
found to depend on a quantity n", which measured the change in angular velocity 
This it was found impossible to evaluate so long as the investigation was confined 
to second-order terms. It now appears that equation (90) is in effect an equation 
determining n n . The equation is linear in n f \ so that it gives only one valine for n ff . 
When n n has this value we are on the true linear series, but if n ,f has any other 
value our solution, when we try to extend it to third-order terms, degenerates into a 
solution of the type of (85), with which no progress can be made. Our plan, then, 
is to evaluate the terms which occur in equation (90) and so obtain the value of n ,r . 
On inserting this into the values of p y q, r, s which were obtained in the previous 
paper, we complete the solution as far as the second-order terms, and can then proceed 
to the stability criterion. 

Numerical Computations. 

17. It is at once apparent that the evaluation of C n , t 12 , . .. ft l5 ft 2 , ... , given by 
equations (58)~(67), can be made to depend on integrals of the same type as occurred 
in the previous paper, namely integrals defined by 



J 



ABC 



>00 



dX 



AABC... 
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In addition to the integrals which were given in the previous paper, # I have 
calculated the following values of additional integrals needed for the computation of 
third-order terms. 



J B5 = 0-9637364, 
J Bao , = 2*948174, 



J 



J A 6 = 

Ja 4 b 2 — 

^A 8 B 8 = 

Ja 8 c s ~~ 

A 2 B 2 C 2 

J 



AB= 



J 



AB'C 3 — 

J B 6 = 

J B S C 5 = 



0-0001716853 : 

0-002877084, 

0*01256355, 

0*0365761, 

01172029, 

0-260331, 

0-802462, 

JL ZiZi J. U J. , 

3-837093, 



A? 



J 



J 



J 



J 



J 47 = 0-0000434660, 

* 

- 0*000757892, 
= 0*00334922, 
= 0*00986287, 
= 0*0317157, 
= 0*0704673, 
= 0*218724, 
= 0*332181, 
: 1*04926, 



A 5 B 2 

J A 4 B 3 
J A 4 C 3 
A 3 B 2 C 2 

Ja 2 b 5 

A 2 B 2 C S 

ti AT16 



AB* 



AB 3 C 3 



J B4C = 1-392786, 
J BC4 = 4*316064, 



J 



A 5 B "™ 

Ja 4 bc = 

J A 3 B 2 C = 

J, 



A 2 B 4 



Ja 2 bc 3 = 

^AB 4 C ^ 
J ABC* = 

J B 5 C = 

J B 2 C 4 = 
J C 6 = 

J A 6 B = 

J Asun = 



A 5 BC 



J 



J 
J 



A 4 B 2 C — 

J A 3 B * = 

A 8 BC 3 = 

A 2 B 4 C = 



J A 2 BC* _ 
J A^sn = 



J 



AB 5 C 

AB 2 C* 
J 



AC S — 



0'0006851396, 

0-004027127, 

0-01783828, 

0-0565289, 

0-1698776, 

0-377266, 

1-177387, 

1781565, 

5*679960, 

12-448855, 

0-000177533, 

0-00106671, 

0-00477536, 

0-0152016, 

0-0460906, 

0-102374, 

0-321581, 

0-485553, 

1-55682, 

3*42068. 



J B80 , = 2-024097, 
J cs = 6-336915, 



J A5n -~ 



A 5 C 



A 4 C» 



"A'BC 8 — 

J A 2 B 3 C = 

T — 

d A 2 C 4 — 

T — 

"AB 3 C 2 ~" 

J AC 5 ~ 

T — 

"B 4 C 2 ~" 

J BC 5 = 



0*000942028 ; 

0"005675881. 

0*02547622, 

0*0811833, 

0*2473229, 

0*549771, 

1*731910, 

2*621427, 

8*391280, 



J A6n — 



A 6 C 

J A 5Q2 

J A 4 BC 2 
J A s B s c 

J A 3 C 4 

J A 2 B 8 C 2 

J A 2Q5 

J AB 4 C 2 
J ABC 5 



0*000245881 

0*00151097, 

0*00684621, 

0*0219023, 

0*0672671, 

0*149566, 

0*473857, 

0*716300, 

2*30256, 



The first task is the evaluation of C 11? C 12 , ... , C 33 . These quantities are connected by 
the three relations (68)-(70), so that only three of the six quantities need have been 
calculated. I have, however, calculated four coefficients, namely, C 12 , C 13 , t 22 and c 33 
quite independently, the object being to obtain a check on the accuracy of the 
computations. In this way a check is obtained at every step of the computations. 

* The following are true values of integrals which were incorrectly given in the previous paper : — 

J BCC = 1-3419631, J AABCC = 0-09510324, I cccc = 0'4781180. 

In each case the error was one of printing only, and did not enter into the computations. The whole 
table given in the previous paper was recomputed before extending it for purposes of the present paper. 
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The values obtained were 

C l2 = -0*00027991L-0*0093206J« + 0*0103775N 

-0*00458151 -0*0016151m +0*004026811 +0*0042388, . . (91) 

C 8 = -0*0003503iL+0*0035301iE-0*0295936N 

-0*00242841 +0*0074994ttt -0*0010296tt +0*0044353, . . (92) 

t 22 = 0'00010103L+0*0126532J1 + 0*0153304N 

-0*02508191 + 0*0010344m -0*0017407tT -0*0015791, . . (93) 

ts8= 0*000136l3L + 0*0062287Jl + 0*035314lN 

-0*02615751 -0*0035230m +0*0007874tt -0*0016798. . . (94) 

Two values of C 23 can now be deduced from equations (69) and (70) respectively. 
These are found to be 

+ 0*0002337 -0*047957 -0*123015 

c 23 = H M W 

+ 0*0002343 -0*047962 -0*123004 

+ 0*16424 -0*001361 -0*001636 -0*003242 

I lit tt 

+ 0*16420 -0*001360 -0*001636 -0*003227. 

The agreement of these values provides a check on the computations of the 
coefficients t 12 , ..., and of the integrals from which they have been calculated. 
In virtue of relations (73), equations (68)-(70) become 

sk+l + ™ = °> W 

a 4 o 4 c 

»£+^-° <-> 

while relations (73) and (74), of which only three are now independent, may be 
represented by 

*»-*^ n =0, ; (98) 

13 2 a b c 4 

c 22 -i--^J»-o (100) 

On substituting for C 12 , C 13 , and C 22 from equations (9l)-(98), these become six linear 
equations for 3L ? JE, N, I, VX 9 and tt. 
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The solution of these equations is found to be 

% = -12*6275, ffl = -0*0307056, N - -0*0044636, 

I = - 0-0116194, m - 0*42602, U = 1*15365. 

These values have been checked by insertion, not only in the six equations from 
which they were directly derived, but also in the remaining equations (73) and (74). 

We now have all the material necessary for the evaluation of % u l£ 2 and lt 3 , 
which, it will be remembered, are the right-hand members of equations (77)-(79) 
respectively. These have been computed independently, and I found 



x = -0*0016803 + 0*228894^', 
a = +0*0026780-0*301110^, 
3 - + 0*0023739-0*385644^". 

These values ought to satisfy (cf. § 16) 

3* 1 + » 2 + *3 = 0, 
in place of which I find 

3% + % + % = 0*000011 -0*00007^, 

but the error is no greater than might reasonably be expected in view of the very 
large number of operations in each computation. # 
The coefficients in equation (90) are found to be 

—k f 2 ~i-Jc 2 - 0*0058753, —.Vi-kh = 0*00024949, 

so that the equation itself becomes 

5*8753 (-0*0016803 + 0*2288940 - 0*24949(0*0026780-0*301110^0, 

and the solution is found to be 

n" - 0*007423. 

Completion of Second Order Solution, 
18. On substituting the value just obtained for n ff into the values for p, g, r, s, 
found in the previous paper (§ 34), I find, 

p = 3*124954, 

q = -0*103164, 

r - -0*015236, 

s = -0*256962, 

thus completing the figure to the second-order terms. 

* Each of the quantities Si, S2, S3 has been computed by expressing it as a sum of integrals of the 
type tabulated on p. 21. The first term in Si, namely 13iKa%~ 4 J AAA , may be thought of as a typical 
term. Each of the quantities Si, S 2 , S 3 consisted of 326 such terms, so that 3Si + S 2 + S 3 is a sum of 978 
such terms, each of which, it must be remembered, is evaluated by a fairly lengthy series of computations. 
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The corresponding rotation is given by 

2 



27T/0 



w + eV = 014200 + 0*007423e J . 



We notice at once that « increases as we pass along the pear-shaped series, whereas 
in Sir G. Darwin's solution w was found to decrease. Thus the present solution 
diverges in essentials from that of Darwin. 

On the other hand the present solution is similar to that for rotating cylinders^ in 
which the rotation was also found to increase as we passed along the pear-shaped 
series, and as we shall now see, the increase is at a very similar rate. 

Our present figure, as far as the first order of small quantities, is 

9 2 o 

x , y z 



+ -~^+-^+e(~0'079x s + 0'l27xy 2 +0'W - 1, 



3*55 0'664 0*424 



while the cylindrical figure, on replacing the parameter lO a/a used in the second half of 
the two-dimensional investigation by e, was found to be 



rp& n »" 



5 0-555 v . y i 

A comparison of these two figures shows that the two e's may be taken to be very 
approximately the same. As regards angular velocity, the value in the present three- 
dimensional problem has been found to be 



w 2 



- = n + n"e* = 0'14200 (l + CT05227e 2 ) ; 



2ltp 

while that in the two-dimensional problems wast 

— = - + 8 ^ - 0*3750 (l + 0*0513e 2 ). 
2tt P 8 448 v ; 

Thus, in so far as it is possible to compare the three-dimensional problem with its 
two-dimensional analogue, we may say that the two rotations are in very close 
agreement. 

Calculation of the Moment of Inertia. 

19. We know that the pear-shaped figure will be stable or unstable according as 
the angular momentum increases or decreases as we pass from the critical Jacobian 
ellipsoid along the series of pear-shaped figures. 

* "On the Equilibrium of Rotating Liquid Cylinders," i Phil. Trans./ A, vol. 200 (1902), p. 67. 

f I have recently repeated the calculations of the two-dimensional problem as far as the third order of 
small quantities, including the evaluation of co 2 ? and find that the numbers given originally for the figure 
and angular velocity are absolutely correct. See, however, § 21 below. 
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The moment of inertia of the pear-shaped figure about the axis of rotation, say 
MF, will be given by 

MF = p (x 2 -\~y 2 ) dxdydz 



where the integral is taken throughout the volume of the pear-shaped figure. We 
have, by our choice of the coefficient s, ensured that the volume of the pear-shaped 
figure shall remain always equal to that of the original ellipsoid, so that we have 



M 



TTp 



abc, 



and therefore 



F 



3 



:7T J 



(x 2 + y 2 ) 



2 \ dxdydz 



abc 



» • * • 



(101) 



Let us transform to co-ordinates x\ y\ z\ given by 



x = 



x 

a 



t __ y 

y - £, 



z 



z 
c 



. . • • . • * 11 U Zi) 



so that the critical Jacobian ellipsoid is reduced to a sphere of unit radius, and the 
pear-shaped figure is reduced to a distorted sphere. With this transformation, 
equation (101) becomes 

jfe a =! fff (ci 2 x f2 + b 2 y f2 ) du f dy r dz r , (103) 



4 



7T J J 



where the integral is taken throughout the figure bounded by the surface 



/ / ™'2 



J2 J2 



x' 2 +y' 2 + z' 2 -l + e^- U~ + /3|^ +yj +k J 



-tie 2 



■Lx'\W\ ,, 

-. — r — tt~ ~r » » * ~t o 



a 



b 4 



(104) 



Let r 2 be written for x/ 2 + y n + z f2 , and let us further put 



x = rx, 



y = ?y> s = ?% 



so that x, y, z, are co-ordinates on a sphere of unit radius. Equation (104) becomes 



X / X 2 . ^y 2 . Z 2 \ , "KiC 



^ 1+ ^a[ a 7 + fi V + y ? +& a 



+ ie 2 



"LxV , My V 4 , 



= (105) 



Let us suppose that r, the radius vector to the boundary of this distorted unit 

sphere is given by 

r = 1 + ef + e 2 g + e 3 h + ... . 
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On substituting this value for r into equation (105) and equating the coefficients 
of e and e 2 , we obtain 



x 



f 2 + 2g + f- 

a 



2f + - 

a 



/ _2 



p o 2 ' 

a 3 ?r c 







(106) 



3 ( x + ^ +y |j + K 



+ 4 



L-« 



x- My 4 



a 



b' 



T 



.4- o 

ess j f .? 



= 0. . (107) 



If c£Q is an element of solid angle, equation (103) may be written in the form 



I 



„2 _ 



3 



e t* 



■IT 



20 



7T J 



(a 2 x 2 + b 2 y 2 ) r^drclQ, 

(a 2 x 2 + by) (1 + 5ef + 5e a g+ 10e*f 3 )dQ. 



Hence we find that F may be written in the form h^ + Ak 2 , where, as far as e 2 , 



F = 



4 



7T J 



(aV+Z>y;(g + 2f 2 )<ZQ. 



The integral is here taken over the sphere of unit radius, and so can be easily 
evaluated. Carrying out the necessary computations, I find 

4/ = i ( a 2 +b 2 ) = 0-844105. 



lX/C 



0'079156e 2 . 



Thus the moment of inertia is given by 

Mk 2 = 0"844105M (l~~0'09378e 2 ) 

20. This again differs from Sir G. Darwin's result, in which it will be remembered 
it was found that the moment of inertia of the pear increased as e 2 increased. But 
we shall now see that the difference between the two results agrees exactly with 
what was to be expected from the different values of n ff used in our two solutions. 

Sir G. Darwin chose his parameter e (which I shall denote by e B ) in such a way 
that the longest radius vector of the pear-shaped figure was a (l + 0'1482e D ). In 
my solution, the longest radius vector is found to be a (l + 0'1309<3j), where e z denotes 
my parameter e. Hence I find as the relation between our parameters 



e 



D 



0'8833e T . 



Darwin's rotation was given by 



00 



2tTTp 



0'14200(l-0-1443066e D 2 )= 014200 (l-0M1259ej 8 ). 



FIGURE OF EQUILIBRIUM OF A ROTATING MASS OF LIQUID. 27 

The general solution, apart from special values for the rotation, is such that, in my 
notation, 



n + eJn" = 014200 ( 1+ — ^- — e 



2tt P v 014200 

Hence I find that Darwin's solution ought to coincide with the solution obtained 
in my previous paper on assigning to n" a value w" D given by 

n"»= -0-015988. 

* 

The solution given in the present paper, which is believed to give a true figure of 
equilibrium, is derived from the general solution by assigning to n" a value n"j given 

by 

n" z = + 0*007423. 

In my previous paper, it was shown that there was, so far as second-order terms 
only were concerned, a doubly infinite series of figures of equilibrium, and it was 
found that these could be defined in terms of two independent parameters, e and £ 
where f was the same thing as e/n". It accordingly appears that any figure of 
Darwin's series differs from my figure of equilibrium having the same value of e 3 
through his f being different from mine. The excess of my f over his will be 

&-d = e?{n" z -n"») - +0'023411e 2 . 

It is readily seen that the increase in the moment of inertia over that of the 
critical Jacobian, which has been called AF, will be a linear function of f and e 2 , say 

/&£-{■ He j , 
Hence the excess of my angular momentum over that of Darwin, will be 

The value of /a is easily obtained by allowing e 2 to vanish in the analysis given on 
p. 72 of my previous paper. The quantity /m then appears as the rate of increase 
of F as we pass along the Jacobian series of ellipsoids. The general value of F, in the 
notation there used, # is 

F = \{al 2 + V 2 ) - ia 2 (l-1271347f)+i6 a (l + 9'20894f), 
whence 

M =^4¥^ = -7*84851. 

The excess of my value of F over that of Darwin ought accordingly to be 

M ^_ D = -01830606/. 

I have changed the sign of £, which, by an oversight, had been taken with the opposite sign in my 
formulae as printed in the previous paper. 

ill M 
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My value of Ah 2 having been seen to be — 0*079 156e/ 5 it follows that Darwin's 

ought to be 

AF = 0-10390e/ - 0'13317e D 2 

so that his value of k 2 as deduced from my calculations ought to be 

F = 0*8441 (l + 0*15777e D 2 ). 

In point of fact the value actually given by Darwin was 

F - 0'8441 (l-f 0*157786e D 2 ). 

It appears, then, that Darwin's moment of momentum agrees exactly with mine, 
as it ought, except for the difference introduced by the different values we have taken 
for n /; . But, besides showing this, the calculations just given provide a check on the 
accuracy of the computations of both of us. Although our figures^ as far as the 
second order, have been calculated by very widely different methods, their moments 
of momentum have been found to agree very closely.^ 

21. The moment of momentum in the cylindrical problem was announced in my 
two-dimensional paper to increase with increasing e 2 . 

On repeating the computations of this paper, I find that the coefficients in the 
equation of the surface were correctly given, but the final computation of F was 
erroneous, t The corrected formula becomes 

F-£ 2 (l-0-lG79e 2 ). 

The Stability Criterion, 
22. We have now found for the pear-shaped figure of equilibrium, 



2 
CO 



= 014200 (l + 0-05227e 2 ) 



2irp 

tf= 0*8441 (l-0*09378e 2 ) 

* It may perhaps be added that before I had discovered this check on my computations, I had calculated 
the mass of my pear-shaped figure by direct integration, using the method of § 19. The total mass ought 
of course to come to exactly M. I find that the terms in a 2 , a/3 . . . result in an increase of mass 
0-0689519M6 2 , while those in L, M, N ... s balance this with a decrease - 0*0689514M^ 2 . 

f The error is in the very last stage of all ; the value of |-y on p. 95 ought to read 

\y = 20-25a 6 -118-93a 4 +41-071a 2 + 229-51, 

and this leads to the formula 

¥ = & 2 (1-167 -88(92) = V(l -0'16788e 2 ). 

I have also calculated ¥ by the method of § 19 of the present paper, and found 

W = & 2 (l-0-16791e 2 ). 



FIGURE OF EQUILIBRIUM OF A ROTATING MASS OF LIQUID. 29 

It follows that the moment of momentum M is given by 

M = M u (l-0 , 06765e 2 ). 

Thus M < M , so that the moment of momentum decreases as we pass along the 
series of pear-shaped figures, and this series is therefore unstable. 

Summary and Discussion of Results. 

23. Throughout the present paper, and my previous paper on the same subject, the 
critical Jacobian ellipsoid which bifurcates into the pear-shaped series of figures of 
equilibrium, has been taken to be 

Any adjacent figure, whether of equilibrium or not, may be supposed to be 

where ^ is a function of x, y, z in which the coefficients are numerically small. For 
special values of <£, this figure will be one of equilibrium. So long as we consider only 
figures which differ infinitesimally from f=0, all the possible figures of equilibrium 
form a linear series, and <j> is of the form 

= eP, (108) 

where P is a function of x, y, z and e is a parameter which must be so small that e 2 

can be neglected. 

In the previous paper it was shown that as soon as e 2 is taken into account, there 

must be supposed to be a doubly infinite series of figures of equilibrium. The general 

form of is 

^eP + e'Q + ^Q 1 (109) 

where £ is a second parameter of the same order of magnitude as e 2 , but capable of 



2 
CO 



varying quite independently of e 2 . The value of - — for this figure of equilibrium is 

greater by £ than the value for the critical Jacobian. The possible figures of 
equilibrium may be thought of as lying inside a rectangle having e, £ as rectangular 
co-ordinates. 

In the present paper I have carried the investigation as far as e 3 , and find that the 
value of <f> as far as third-order terms must be of the form 

= eP + e 2 Q + ^Q 1 + e 8 (R + KP) . ...... (110) 

where R is a new function of x. y. z and K is a constant. At first sight K appears to 
be at our disposal, for if we replace the parameter e by a new parameter e + 0e 3 , we can 
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vary K as much as we please. But on examining the problem in detail it is found that 

K is always infinite except for one special value of £ For this special value of f we 

can eliminate K altogether by a new choice of parameter. But for other values of £ our 

solution is only valid if e 3 K is small, and, replacing e 3 K by a new parameter e 1 , the 

solution reduces to 

<p = O.P+O.Q + ^P 

and so returns to the original first-order solution (108). 

Thus, except for one special value of £ it is impossible to carry the linear series 
outside the second order rectangle ; if we attempt to do so, the solution keeps lapsing 
back into the first-order solution. 

In the previous paper I gave reasons for believing that Sir G. Darwin, in his well- 
known investigation of this problem, had introduced a spurious equation of equilibrium. 
This extra equation could only be satisfied by assigning to £ a special value, namely 

f = -0'015988e 2 . 

Sir G. Darwin accordingly gave this value to £ so that the value of w decreased on 
passing along his series of pear-shaped figures, and, assuming this value for £ he showed 
his series to be stable. 

But the investigation of the previous paper showed that there was no need to assign 
this special value to £ and the present investigation has further shown that with this 
value of f it is impossible to extend the series beyond second-order terms at all. 
There is only one value of f which leads to a real linear series of configurations, and 
this is shown in the present paper to be 

£^+0'007423e 2 . 

Thus as we pass along the true linear series o? continually increases. The angular 
momentum is however found to decrease, so that the pear-shaped figure is shown to be 
unstable. 

24. The amount of computation involved in the problem has proved to be very 
great, and as the whole question of stability or instability depends on the sign of a 
single term at the end of all this computation, the question of numerical accuracy 
becomes one of great importance. 

The difference between my second-order figure and that of Darwin arises solely 
from the difference in the value of £ The moment of inertia of such figures is a 
linear function of £ and a very simple calculation gives the rate at which it ought to 
vary with £ Allowing for this difference in £ I find that my computations give for 
the moment of inertia of Darwin's figure (in terms of Darwin's parameter e D ), 

MV(l + 015777e D a ), 

while Darwin calculated as the value of the same quantity 

MV(l + 0'157786e D 2 ). 
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Since two independent sets of computations, conducted by entirely different methods, 
have been found to lead to precisely the same result, it seems highly probable that 
this result is accurate. The agreement just mentioned may reasonably be regarded as 
guaranteeing the accuracy of all the second-order computations, both of Darwin and 
myself. 

The actual criterion of stability, however, depends on the value assigned to f, and 
since this depends in turn on the third-order terms, no check by comparison with 
Darwin's work is possible. Some support is given to my value of f by comparison 
with a parallel investigation of the "Equilibrium of Eotating Cylinders," which I 
published some years ago. Adjusting the parameters so that e shall have, as closely 
as possible, the same physical interpretation in the two problems, I find for the factor 
expressing the increase of co 2 as we pass along the series of pear-shaped figures : 

1-f 0*05227<3 2 for the three-dimensional problem, 
1-f 0*051 3<3 2 for the two-dimensional problem. 

Apart from this, the checks I have used in the present paper are such that I believe 
it would have been very difficult for any error to escape detection. 

25. The main object of the paper is achieved as soon as the pear-shaped figure is 
shown to be unstable. It is at the same time of interest fco examine the bearings of 
this result on the wider question of which it is a part. 
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In fig. 1 let PP 7 represent the series of Jacobian ellipsoids, the part PO (drawn 
thick) representing the stable part of the series, and the part OP 7 (drawn dotted) 
representing the unstable part, so that is the point of bifurcation. Let a diagram 
be drawn about this line having the angular momentum always represented by the 
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vertical co-ordinates, so that for instance all systems, whether in equilibrium or not, 
which have the same angular momentum as the critical Jacobian ellipsoid, must be 
represented in the horizontal plane through (this plane having as many 
" dimensions " as are necessary). In this diagram the pear-shaped series of figures 
will lie below this plane in the neighbourhood of 0. 

There are two a priori alternatives, represented on the right-hand and left-hand sides 
of fig. 1. In the first place, it is possible that after passing a certain distance OB/ 
along the pear-shaped series with decreasing angular momentum, we come to a region 
R'S' in which the angular momentum again increases. Any horizontal line in the 
diagram ought, on the principle of stable and unstable configurations of equilibrium 
occurring alternately, to meet stable and unstable branches of the linear series 
alternately. Thus the branch B/S' ought to be stable, so that R' would be a 
point of bifurcation on this series, and the other series through R', say 1VT, would 
be unstable. 

The alternative possibility is that the pear-shaped series of figures proceeds 
continually downwards in the diagram, so that its angular momentum continually 
decreases. 

26. Either of these two possibilities removes a difficulty to which Sir G. Darwin 
has drawn attention.^" In what is commonly referred to as Roche's problem an 
attempt is made to discover the other end of the pear-shaped series of figures, on the 
supposition that this other end represents two detached masses revolving about one 
another. Two such series of figures appear to existf ; in one the satellite is spheroidal 
except for tidal distortion, while in the other it is ellipsoidal. The former series has 
been shown to be stable, the latter unstable. 

As the angular momentum decreases on passing along these series, the distance 
between the two masses also decreases until a point is reached at which the two 
series coalesce, the configuration of bifurcation being one in which the radius vector 
from the centre of the primary to that of the satellite is equal to 2'457 radii of the 
primary. J If the distance between the masses is decreased still further, the remaining 
configurations form an unstable series. Sir G. Darwin found a difficulty in the 
instability of this series, since he believed it to be the far end of the pear-shaped 
series which he thought stable. We now see, -however, that this series may, without 
change of stability, join up with either the series TR' on the right-hand of our 
diagram, or the series RO on the left. 

* 'Coll. Works,' III., pp. 515-524. 

t Roche's problem has only been solved strictly by imposing sphericity on the primary and assuming 
the satellite to be infinitesimal. Sir G. Darwin's work (* Coll. Works,' III., p. 436), leaves little room 
for doubt that Koche's result may be extended in the way I have stated. 

\ If the satellite is not infinitesimal, the radius vector depends on the ratio of the masses, but always 
lies between the narrow limits 2*457 and 2*514 times the radius of the primary (see Darwin, loc. cit., 
p. 507). 
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27. We must now consider what motion is to be expected in a Jacobian ellipsoid 
which has reached the point of bifurcation at which instability sets in. Poincake 
remarks # that if the pear-shaped figure proved to be unstable, " la masse fluide 
devrait le dissoudre par un cataclysme subit." 

After reaching the point in our diagram, the mass cannot move along the 
pear-shaped series, since this would involve a decrease of angular momentum. It 
may be thought of as moving along the unstable branch OP' of the series of Jacobian 
ellipsoids for an infinitesimal time until some slight disturbance brings its instability 
into play. - 

Now of all the vibrations of this figure, it is known that one only is unstable, 
namely that corresponding to the third zonal harmonic of the ellipsoid. The initial 
motion of the figure must then be one in which the displacement at every point of 
the surface is proportional to the third zonal harmonic. 

Thus the fluid begins by describing exactly the pear-shaped series, but as soon as 
the changes in angular momentum become appreciable, it leaves this series, and passes 
through a series of configurations represented in the region above in fig. 1 . These 
may at first be thought of as lying parallel to the pear-shaped series, but above O. 

If there is a stable branch such as B/S ; which ultimately passes above 0, it is 
conceivable that the series of non- equilibrium configurations might ultimately coalesce 
with the series of equilibrium configurations B/S', and the motion would be continued 
along this series. In this case, M. Poincaee's " cataclysme subit " would consist in a 
jump from the stable series PO to the stable series R/S ; . 

Judging from the results of my parallel investigation on the configurations of 
rotating cylinders, this possibility does not seem at all likely. It is, I think, much 
more probable that the pear-shaped series lies like the series OR in my figure. 

In this case also the liquid would move through a series of configurations which 
would initially be close to the series OR, but would get continually father removed 
from configurations of equilibrium. The protuberance resulting from the initial third 
harmonic displacement would develop in a manner somewhat similar to that of the 
pear-shaped figure, but as the motion would necessarily be possessed of a considerable 
amount of kinetic energy, the phenomenon would be a dynamical and not a statical 
one. If the configuration represented at in fig. 1 is the highest stable configuration 
possible for a single mass of liquid, this kinetic process can end in only one way, 
namely in the separation of the mass of liquid into two parts. As the third harmonic 
displacement develops, the region of the pear which moves with greatest velocity is 
known to be the extreme end of the protuberance. It is, therefore, natural to suppose 
that this part of the figure would be shot away first. Moreover as the departure 
from a figure of equilibrium is probably pretty pronounced before the separation takes 
place, it is likely that the mass in question will be shot away with a considerable 
velocity. 

* Letter to Sir G. Darwin in the latter's ' Coll. Works/ III., p. 315. 
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After this projectile has left the main mass, the angular momentum of the latter 
(measured of course per unit mass) will be reduced, and as the tidal influence of the 
newly-born satellite is gradually withdrawn, the primary may settle down to a state 
of stable equilibrium in which its figure is again that of a Jacobian ellipsoid. A 
cycle of processes such as this might very conceivably constitute the life- history of a 
rotating body after once it had passed the critical state represented by the point of 
bifurcation on the Jacobian series. There seems to be no reason why the protuberance 
should always develop at, and be shot off from, the same end of the Jacobian 
ellipsoid. 

Thus it appears that the instability of the pear-shaped figure leads us to contemplate 
a series of processes of much the same nature as would have been expected if the pear 
had proved to be stable, except that we are now led to assign a much shorter time to 
these processes. The problem has ceased to be one of statics and has become one of 
dynamics ; the phenomenon is no longer one of slow secular change but of collapse 
and explosion. The mechanism of planetary birth which is now suggested is so rapid 
that there need be no difficulty in supposing a planet to have been through the cycle 
several times ; had the pear proved to be stable the cycle would probably have been so 
slow as to create a real difficulty. 



